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ABSTRACT 

Two coalescing black holes (BHs) represent a conspicuous source of gravitational waves (GWs). The merger 
involves 17 parameters in the general case of Kerr BHs, so that a successful identification and parameter ex- 
traction of the information encoded in the waves will provide us with a detailed description of the physics of 
BHs. A search based on matched-filtering for characterization and parameter extraction requires the develop- 
ment of some 10 15 waveforms. If a third additional BH perturbed the system, the waveforms would not be 
applicable, and we would need to increase the number of templates required for a valid detection. In this letter, 
we calculate the probability that more than two BHs interact in the regime of strong relativity in a dense stellar 
cluster. We determine the physical properties necessary in a stellar system for three black holes to have a close 
encounter in this regime and also for an existing binary of two BHs to have a strong interaction with a third 
hole. In both cases the event rate is negligible. While dense stellar systems such as galactic nuclei, globular 
clusters and nuclear stellar clusters are the breeding grounds for the sources of gravitational waves that ground- 
based detectors like Advanced LIGO and Advanced VIRGO will be exploring, the analysis of the waveforms 
in full general relativity needs only to evaluate the two-body problem. This reduces the number of templates of 
waveforms to create by orders of magnitude. 



1. INTRODUCTION 

The detection of merging black holes (BHs) is the holy grail 
of ground-based detectors of gravitational waves (GWs) such 
as LIGO and VIRGO. For the search for GWs of compact bi- 
naries, the availability of accurate waveform models for the 
full merger is crucial. Thanks to the success of numerical rel- 
ativity in simulating the l ate inspiral, me r ger and ringdown 
of a binary of two BH s (Pretorius 2005; Campanell i et al.| 
2006 Baker et al. 2006| l, we are now able to perform a search 
with realistic templates. The conjunction of post-Newtonian 
modelling of the inspiral phase and full numerical relativis- 
tic simulations of the merger and ringdown is now a reality 
for comparable- mass scenarios of mass ratio s up to ab out 10 
(see e. g. Buona nno & Da mour 1999; Buonanno et al.||200"7*t 
Ajith et al.||2009| |San tamarfa et al.||2010[ l Nevertheless, the 
high cost or full numerical relativity simulations constitutes a 
serious limitation to the development of waveforms. 

It has been estimated that we need a bank of 10 13 wave- 
forms for the identification of a binary of two B Hs with an 
F-statistic b ased grid search at a minimal match ( |Cornish ~&\ 
|Porter||2~007| ). On the other hand, since the BHs are most 
likely Kerr, the number goes up to ~ 10 15 to cover the spin 
parameter. In the case of stellar-mass BHs, since we can- 
not detect the sky-location of the source, we are limited to 
a two-dimensional space, the masses of the BHs, so that the 
number is considerably reduced; we need to create about 
~ 10 5 - 10 6 templates. This is the reason why in the last 
years there has been a significant effort in developing alter- 
native approaches, such as Monte-Carlo schemes, genetic al- 
gorithms, Metropolis-Hastings methods and Nested Sampling 
techniques ( Cornish & Porter]|2006 Lang & Hughes|[2006| 



|Gair & Porter|2009||Petiteau et al.|2 009 ). 
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Campanelli et al. (2008 ) addressed for the fist time fully rel- 
ativistic long-term numerical evolutions of three equal-mass 
BHs and found that the merger dynamics is very distinct from 
binaries. In particular, the trajectories were intricate and led 
to singular waveforms, as e.g. their figure 4 shows, in which 
we can see two mergers. Recently there has been an effort 
in calculating in detail the wave forms of systems of three and 
four BHs interacting in full GR. |Galaviz et al.| ( |2010| l have de- 
veloped a knowledgeable scheme to study the waveforms of 
such configurati ons and find int ricate templates for the waves. 
Also, |Jaramillo & Lo usto (2010) have addressed the problem 
of critical BH separations for the formation of a common ap- 
parent horizon. The authors study in detail the aligned equal 
mass cases for up to 5 BHs. 

If we increase the number of BHs involved in the GW, the 
number of templates to develop increases enormously. Putting 
it in Neil Cornish' words, "The sensitive dependence on ini- 
tial conditions will send the template count through the roof". 
It is consequently important to understand the limits imposed 
by the physical systems which harbour these sources of GWs. 
Therefore we address the question of the existence of a system 
with more than two BHs in a relativistic regime. In section [2] 
we calculate the probability of having a relativistic three-body 
encounter in a dense stellar cluster with three BHs initially 
unbound. In section [3] we estimate the possibility that an al- 
ready formed binary of two BHs interacts relativistically with 
a third BH. In section|4]we summarise our results and give the 
conclusions. 

2. THREE-BODY RELATIVISTIC ENCOUNTERS OF UNBOUND BHS 

Rough estimtes will be sufficient to show how unlikely 
triple relativistic encounters are. Therefore, for simplicity, 
we assume that all BHs in a given stellar system have the 
same mass, m. Let us assume that we have two of them fly 
by with a periapsis distance of a few Schwarzschild radii, 
^peri = dGm/c 2 . Therefore, d is the periapsis distance in units 
of Gm/c 2 , with G the gravitational constant and c the speed of 
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light. If the relative velocity between the two objects at large 
separation is V K \, the cross-section for an encounter with pe- 
riapsis distance within d can be estimated as 
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We have assumed that the BHs are in a non-relativistic stellar 
cluster, as observed, V re \ « a <c. Here, a is the 3D velocity 
dispersion in the cluster, defined such that the total kinectic 
energy in the cluster is jMer 2 , where M is the total mass of 
the cluster. 

Hence, the encounter rate for one BH is l/f en c = nSV le \ = 
2irdG 2 m 2 n/(c 2 (j), where n is the number density of BHs. 
In order to obtain the rate of relativistic three-body encoun- 
ters, one multiplies by the probability of having a t hird object 
within a volume 
|Hut|2003D 
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stars during hyperbolic 2-body encounters. Relaxation drives 
the over all expansion and eva poration of the cluster (see for 
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ation time is approximately ( |Binney & Tremaine] 2008 ) 
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Hence, we can estimate the total number of relativistic 3-body 
encounters over the lifetime of a self-gravitating cluster, 
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This is the rate per object. We now calculate the rate for a 
whole cluster of N BHs, 
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with M = Nm and n^iN/R 3 . We ascertain now that the cluster 
is self-gravitating, so that GM/R = a 2 , where R is a typical 
length-scale of the cluster, for instance its half-mass radius, 
i.e. the r adius of sphere containg o ne half of the BHs (see for 
instance [Binney & Tremaine]2 008 ) 
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This equation makes it clear that the probability of even just 
one triple encounter in a relativistic regime (small value of d) 
is extremely low, even if one manages to pack a million stellar 
BHs within a sphere with a radius of 0.1 pc. 

So far, we have considered only a self-gravitating cluster 
of compact objects (COs), such as the stellar BHs, that have 
accumulated at th e centre of a globular cluster as a result of 
mass segregation (Baumgardt et al. 2003a b 2004). However, 
a more promising environment for the accumulation of a large 
number of COs in a small volume is the centre of a galac- 
tic nucleus. There the gravitational force is likely dominated 
by a massive black hole with a mass M M bh ^ 10 s M Q (|Ko- 
rmendy & Richstone | |l995| |Magorrian et al.||1998| [Gultekin 
et al. 2009 ). In that case, a = \/GMmbh/R, and equation |4 



must be modified with a factor \/M/Mmbh- This shows that 
relativistic triple interactions are even less probable within 
COs orbiting MBHs than in a self-gravitating cluster. 

The lifetime of an isolated self-gravitating cluster is lim- 
ited by 2-body relaxation, i.e. the exchange of energy between 



3. THREE-BODY INTERACTIONS BETWEEN A BINARY AND A 
SINGLE BH 

The possibility of a relativistic encounter between three 
single objects being ruled out, we turn to the possibility of 
achieving a 3-body relativistic interaction through the en- 
counter between a binary and a single object. Here the dif- 
ficulty is that the binary itself needs to be already in a rela- 
tivistic regime. This implies that it is emitting gravitational 
radiation at a high rate, hence his lifetime is necessarily very 
limited. Consider a binary of two BHs of semi-major axis a 
and masses m\ and 1112 and a single BH of mass mj, passing at 
a distance of d from the center of mass of the binary. 

Approximately, the binary will need a time few to shrink its 
orbit due to the emission of gravitational waves as estimated 
by |Peters| ( (T964) , 
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Let us now estimate the timescale t-} for the system to in- 
teract with a third black hole that flies by at a distance of < d 
of the centre-of-mass of the binary. If the mass of the third 
BH is ;«3 and the relative velocity to the binary V le i, the cross 
section taking gravitational focusing into account is 
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Therefore, the timescale is 
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We assume that the three objects have a similar mass m; i.e. 
mi « m2 ~ W3 := m. Hence, 
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5 Energetically, the overall cluster expansion can occur only because a very 
small number of stars, at the centre of the cluster, are getting more and more 
bound together. A single tight binary can drive this expansion by releasing 
energy to passing-by objects until it is ejected from the cluster or merges. In 
the next section, we consider whether relativistic encounters can occur during 
interaction between a binary and a single object. 
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We again measure distances in units of GM/c 2 , with a := 
a/(Gm/c 2 ) and d := d/(Gm/c 2 ). Typical interesting values 
are a « a few and d w few tens. Therefore, 
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For a self-gravitating cluster, we can rewrite this in terms of 
N, and R, using V re i w a, 
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This quantity can be interpreted as the probability that a rel- 
ativistic binary has a close encounter with a single object be- 
fore it merges. The total number of mergers that can occur 
in the evolution of a cluster of N objects is A^ merg < N, with 
A'merg = N only possible if all objects merge together (a sce- 



nario toforma MB H from a cluster of stellar BHs, see |Lee| 
1993 Kupi etal.|2006| l. Hence, the total number of relativistic 
single-binary interactions over the lifetime of a cluster (evolv- 
ing through a succession of binary mergers) is 
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4. DISCUSSION 

In this work we have addressed the formation of systems 
of three BHs in the strong gravity regime. For that we have 
first studied the probability that three BHs interact in a dense 
stellar cluster and we conclude that it is totally negligible. We 
have then addressed the possibility that a binary of BHs which 
has previously formed in the cluster interacts relativistically 
with a third BH. We judge that the stellar system harbouring 
the BHs needs to have unachievable densities. 

Based on simple physical arguments, we have established 
that the time scales for a triple relativistic encounter to occur 
in a cluster of stellar-mass compact objects is extremely long. 



We recall that we consider two extreme cases. In the first sit- 
uation, we neglect the existence (and formation) of binaries. 
Hence, we consider that three single objects have to find them- 
selves, by chance, within a few (tens of) Schwarzschild radii 
of each other. In this case, we have admitted that the clus- 
ter cannot survive (with a high stellar density) for more than 
about 100 relaxation time. Indeed, in such a long time scale, 
most of the cluster would expand to lower and lower densities 
and a very significant fraction of the object would escape, with 
only a very small number of object s getting closer and closer 
to provide energy for this evolution (Heggie & Hut 2003]f. The 
second (much idealized) case is that of a cluster made mostly 
of binaries. In that case, one can hope that it would suffice 
for an object to interact closely with a binary but in order for 
the interaction to be relativistic for the three objects, the bi- 
nary must be so tight that its lifetime is limited by emission 
of gravitational waves. Accordingly, we consider that the life- 
time of the whole cluster is limited by the successive merger 
of binaries. Any real cluster would present a situation which 
is somewhat in between these two extremes. In particular, the 
evolution of a cluster made of single objects would naturally 
lead to the formation of binaries during core collapse (see e.g. 
Heggi e~& Hut||20"03"] ). We stress that we have made a large 
number of simplifications in our estimates but always in such 
a way as to overestimate the rate of triple relativistic events. 
For instance, we have assumed that all the binaries in a cluster 
are relativistic. This limits their lifetime but non-relativistic 
binaries are useless for our purpose. 

With figure [T] we can estimate what conditions are required 
for at least one such encounter to take place during the lifetime 
of the cluster. For this figure, we have assumed m = 10M Q 

and d = d = 10 4 . The latter values correspond to encoun- 
ters that are only weakly relativistic. Even with such val- 
ues, the figure shows that most of the parameter space for 
single ^ 1 or A^ bin ^ 1 is excluded, either because the clus- 
ter, as a whole, would be a massive black hole (for large N 
and small R) or because the cluster would have such a small 
2-body relaxation time that it wouldn't have time to form (for 
smaller N and small R). Indeed, the evolution of massive stars 
into compact objects (neutron stars or black holes) requires at 
least 3 million years . Because of these constraints, it appears 
that clusters hosting triple relativistic encounters have to be 
made of at least 10 8 compact objects concentrated within a 
region smaller than 10~ 3 pc. In fact, for N ~ 10 s , the size has 
to be of order 10~ 4 pc. The corresponding number density is 
comprised between 10 17 and 10 19 pc" 3 . 

Such values are beyond observed ones by many orders 
of magnitude. For instance, the stellar density in Galactic 
globular clu sters is, at most, of the order of 10 5 pc" 3 (see 
|Harris||1996|and the 2003 on-line update http://www. 
Iphysics . mcmaster . ca/~harris/mwgc . dat) . Al- 
ready the necessary number of compact objects is much larger 
than what one can expect in the kind of clusters that exist. 
A globular cluster might contain up to 10 7 stars but only a 
very small fraction of them would become BHs. This num- 
ber fraction / depends on the initial mass function (IMF) of 
high m ass stars. Some galactic nuclei seem to be top-heavy 
(Maness et al. 2007), so that in principle / ranges between 
10~ 3 and 10~ z , yielding, at most, N = 10 5 BHs in very large 
globular cluster. 

A galaxy contains a much larger number of compact objects 
but most of them inhabit regions of very low density. BHs 
born in the central regions are likely to accumulate at the cen- 
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FIG. I. — Parameter space for the occurrence of triple 
relativistic encounters. For this figure, we consider a self- 
gravitating cluster made of a given number of identical com- 
pact object with an individual mass of 1OM0. The line la- 
beled bin ~ 1 is based on equation 15 with d = a = 10 4 . 
Below this line, in our simple estimate, at least one binary- 
single encounter in which the three objects are relativistic 
should occur during the life-time of the cluster (which is 
limited by successive mergers, see text). The line labeled 
^3, single ~ 1 is based on equation 6lwith d = 10 4 . Below this 
line, one can expect at least one triple relativistic encounter 
between single objects to occur during the life-time of the 
cluster (limited by 2-body relaxational processes). The large 
gray area in the bottom-right section is excluded from the 
parameter space because the whole cluster would be smaller 
than its Schwarzschild radius, i.e. it would immediately col- 
lapse into a single massive black hole. The other gray region 
correspond to clusters in which the 2-body relaxation time 
would be shorter than 10 4 years. Such a cluster would proba- 
bly not last long enough for stars to form and turn into com- 
pact objects. The dashed diagonal lines in gray indicate an 
estimate of the number density in the cluster, from the rela- 
tion n N/R 3 



tre, in the galactic nucleus, through the process of mass segre- 
gation, which is basically an effect of dynamical friction. But, 
in the case of our Galaxy, dynamical friction is only effective, 
over a Hubble time, within a few parsecs of the centre. Hence, 
at most of or der 10 4 BHs may have gathered within the inner - 
most 0.3 pc ( Freitag et al.||2006 Amaro-Seoane et aL)|2010| i. 
A dark mass concentration weighing 4 x 10 b M has been de- 
tected at the Galactic centre through the analysis of the orbits 
of bright IR stars, the so- called S- or SO-stars, around the 
weak radio source Sgr A* (Eisenha uer et al.|2005||Ghez et al. 



[2005] [2008) |Gillessen et al.||20(») |Genzel et al.pOlO) . It is 
generally assumed that it is a massive black hole although the 
only strict constraint on its size is that it has to fit within the 
periapse of the IR stars' orbits, the tightest of which is that of 
S-2 (SO-2), imposing R<5x 10~ 4 pc. If we assume that this 
object is actually a cluster of 10M o BHs, with N sa 4 x 10 s , it 
would need to be so compact in order to host triple relativistic 



encounters that its relaxation time would be extremely short, 
making its existence at the present time an e xtraordinary co- 
incidence (See also Maoz 1998) |Miller|2006| ). Furthermore, 
no mechanism to form such a dense cluster is known. 

While we have limited the analysis to a system of three 
BHs, it is nevertheless obvious that for more BHs the event 
rates are much more unlikely. We therefore conclude that the 
waveforms to develop need only to include two BHs for the 
searches of GWs in the data streams. 
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